Organization of ecosystems in the vicinity of a novel phase transition 
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It is shown that an ecosystem in equilibrium is generally organized in a state which is poised in 
the vicinity ol a novel phase transition. 
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An ecological community consists of individuals of dif- 
ferent species occupying a confined territory and sharing 
its resources [llSS 13- One may draw parallels between 
such a community and a physical system consisting of 
particles. In this letter, we show that an ecosystem can 
be mapped into an unconventional statistical ensemble 
and is quite generally tuned in the vicinity of a phase 
transition where biodiversity and the use of resources are 
optimized. 

Consider an ecological community, represented by indi- 
viduals of different species occupying a confined territory 
and sharing its common resources, whose sources (for ex- 
ample, solar energy and freshwater supplies) depend on 
the area of the territory, its geography, climate and en- 
vironmental conditions. The amount of these resources 
and their availability to the community may change as a 
result of human activity and natural cataclysms such as 
climate change due to global warming, oil spills, defor- 
estation due to logging or volcanic eruptions. 

Generally, species differ from each other in the amount 
and type of resources that they need in order to survive 
and successfully breed. One may ascribe a positive char- 
acteristic energy intake per individual, Sk, of the A:-th 
species (0 < eo < £1 < ■ • ■ )i which ought to depend on 
the typical size of the individualja, la, S lal- The indi- 
viduals of the A:-th species play the role of particles in a 
physical system housed in the energy level Sk- 

We make the simplifying assumption that there is no 
direct interaction between the individuals in the ecosys- 
tem. Our analysis does not take into account predator- 
prey interactions but rather focuses on the competition 
between species for the same kind of resources. 

In ecology, unlike in physical systems where one has a 
fixed (average) number of particles and an associated av- 
erage energy of the system, one needs to define a new sta- 
tistical ensemble. One may define the maximum amount 
of resources available to the ecological community to be 
Emax- The actual energy used by the ecosystem, (E), 
can be no more than Emax 
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smaller the T, the larger the amount of energy utilized by 
the system. Thus, in order to make best use of the avail- 
able resources, a system seeks to minimize this imbalance. 
Second, T behaves like a traditional 'temperature' in the 
standard ensembles in statistical mechanics 9] and, in an 
ecosystem, controls the relative species abundance. 



I. SKETCH OF THE DERIVATION 

Consider the joint probability that the first species has 
rii individuals, the second species has 712 individuals and 
so on: 

s 
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where S is the total number of species (we assume that 
5^1) and 8(u) is a Heaviside step-function (defined to 
be zero for negative argument u and 1 otherwise) which 
ensures that the constraint is not violated. Here, P(n) 
represents the probability that a species has n individuals 
in the absence of any constraint and is the same for all 
species. Note that P{n) does not have to be normalizable 
- one can have an infinite population in the absence of 
the energy constraint. 

As in statistical mechanics ,9i|, the partition function 
(the inverse of the normalization factor), Q, is obtained 
by summing over all possible microstates (the abun- 
dances of the species in our case): 
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On substituting Eq. |(5J into Eq. lO and representing the 
step-function in the integral form one obtains 
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and, strikingly, as shown later, the energy imbalance, 
T ^ 0, plays a dual role. First, for a given Emax, the 



dz 

2'Kiz 



^E^ax-Ylk h{zek) 



(4) 



where e~^'^^^ = ^„ e~^^P{n) and the contour 7 is paral- 
lel to the imaginary axis with all its points having a fixed 
real part zq (i.e. 2; e 7 <^ z = zq + iy, —00 < y < +00 
). The integral is independent of zq provided zq is 
positive J£|. 

We evaluate the integral in Eq. Q by the saddle point 
method'lO'l by choosing zq in such a way that the max- 
imum of the integrand of Eq. Q) occurs when y = 0: 
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where the prime indicates a first derivative with re- 
spect to the argument. Note that the rhs is sim- 
P^y X]fc^fc("fe) with the average taken with the weight 
P(nfc)exp[-zo£fcnfc]. 

Comparing Eq. jSJ with Eq. (^|, one can make the 
identification zq = 1/T and therefore 
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This confirms the role played by the energy imbalance as 
the temperature of the ecosystem. Indeed, the familiar 
Boltzmann factor is obtained independent of the form of 
P{nk). Note that Eq. © with P{n) = 1 leads to a system 
of non interacting identical and indistinguishable Bosons, 
whereas P{n) = l/nl describes distinguishable particles 
obeying Boltzmann statistics. It is important to note 
that terms neglected in the evaluation of the integral in 
Eq. (01 contribute for finite size systems but these vanish 
in the thermodynamic limit. 

In order to derive an expression P{n), consider the dy- 
namical rules of birth, death and speciation in an ecosys- 
tem. In the simplest scenario [ll|, the birth and death 
rates per individual may be taken to be independent of 
the population of the species, with the ratio of these rates 
denoted by x. Furthermore, when a species has zero pop- 
ulation, we ascribe a non-zero probability of creating an 
individual of the species (speciation). Without loss of 
generality, we choose this probability to be equal to the 
per capita birth rate. 

One n iay write down a master equation for the dyna- 
mics [13, U^ uM ^^d show that the steady state proba- 
bility of having n individuals in a given species, P{n), is 
given by the distribution: 



P(n)=P(0) — ,n = l,2,3,... 
n 



(7) 



When X is less than 1, this leads to the classic Fisher log- 
series distribution ; 1 5j for the average number of species 
having a population n, 0(?t.) oc P{n). 

On substituting Eq. Q into Eq. ©, one obtains 
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where the term [a;exp(— efe/T)]"*/nfc is replaced by 1 
when nk = 0. Note that this leads to an effective birth 



to death rate ratio equal to a;exp(— Sfc/T) < 1 for the 
fc-th species. In a non-equilibrium situation, such as an 
island with abundant resources and no inhabitants, the 
ratio of births to deaths can be bigger than 1 leading to a 
build-up of the population. In steady state, however, the 
deaths are balanced by births and speciation (creation of 
individuals of new species). 

It is interesting to consider an ecosystem with an addi- 
tional ceiling on the total number of individuals, Nmax, 
that the territory can hold. In analogy with physics, one 
may define a chemical potential^], /i ^ 0, so that its 
absolute value is the basic energy cost for introducing 
an individual into the ecosystem. Thus the total energy 
cost for introducing an additional individual of the A:-th 
species into the ecosystem is equal to £fc — ^ - effectively, 
all the energy levels are shifted up by a constant amount 
equal to this basic cost. 

The chemical potential may also be defined as the neg- 
ative of the ratio of the energy imbalance to the popula- 
tion imbalance: 
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where (A^) is the average population. It follows then that 
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where a = exp{fi/T) ^ 1. This equation has the same 
structure as Eq. (^. Interestingly, the link between the 
population imbalance and the chemical potential can also 
be established formally starting from Eq. Q, but with 
an additional ceiling on the total number of individuals. 
The introduction of the ceiling on the population leads to 
an additional suppression of the effective birth to death 
rate ratio which now becomes aa;exp(— Efe/T). 

Following the standard methods in statistical me- 
chanics 9J , one can straightforwardly deduce the ther- 
modynamic properties by taking suitable derivatives of 
F = — TlnQ, the free energy: 

F = -T^ ln[l - ln(l - aa;exp(-efc/r))]. (11) 



The average number of individuals in the fc-th species. 
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, is given by 
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In ecological systems, one would expect, in the sim- 
plest scenario, that there ought to be a co-existence 
of all species in our model with an infinite population 
of each when there are no constraints whatsoever or 
equivalently when E^ax = ^max — 00. Noting that 
a = 1 when Nmax = cxd, this is realized only when 
1 — a;exp(— efe/r) = for any k, which, in turn, is valid 



if and o nly ii T = oo and x = 1. We will restrict our 
analysis [l^ in what follows to the case of x = 1. 

Following the treatment in physicsQ, we postulate 
that the number of energy levels, or equivalently the 
number of species, per unit energy interval (the density 
of states) is proportional to the area of the ecosystem 
and additionally scales as e'^ with d > in the limit of 
small e. This is entirely plausible0,|a,0l3 because one 
would generally expect the density of states to have zero 
weight both below the smallest energy intake and above 
the largest intake and a maximum value somewhere in 
between. 

In a continuum formulation, one obtains the following 
expressions for the average energy {E) and population 
(N) of the ecosystem: 
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Note that, except for the second factor in the denomi- 
nator, which is subdominant, these integrals {Io{a) and 
/i(a)) are identical to those of a Bose system. The key 
point is that (in a Bose system and here) they are both 
when a = and monotonically increase to their separate 
finite maximum values at a = 1. 



II. RESULTS AND CONCLUSIONS 

For a given E„iax = Em, Eq. ^ can only be satisfied 
over a finite range of temperatures. The temperature 
cannot exceed Tmax = ^j\/, because (E) cannot be neg- 
ative. At this temperature, (E) — a ^ (N) — Nmax = 
and the territory is bereft of life. Also, the lowest attain- 
able temperature, Tmm, corresponds to Nmax = oo and 
satisfies the equation 



Em - Tmin + T'mt?^l(l) 
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(recall that /i is largest when a — \). 

The increase of {E) on decreasing T is counter-intuitive 
from a conventional physics point of view. The simple 
reason for this in an ecosystem is that a decreasing T 
corresponds to a decreasing energy imbalance (first term 
on rhs of Eq. iQJ thereby leading to a corresponding in- 
crease in the second term, which is {E). This increased 
energy utilization, in turn, leads to an increase in the 
population of the community. 



FIG. 1: Comparison of the results of computer simulations of 
an ecosystem with theory. We consider a system with 100, 000 
energy levels with Ck ~ fc^", k — 1..100, 000, corresponding to 
d — 1/2. We work with a constant Emax (the figure shows the 
results for several values of Emax) and consider a dynamical 
process of birth and death. We have verified that the equi- 
librium distribution is independent of the initial condition. 
At any given time step, we make a list of all the individu- 
als and the empty energy levels. One of the entries from the 
list is randomly picked for possible action with a probability 
proportional to the total number of entries in the list. Were 
an individual to be picked, it is killed with 50% probability 
or reproduced (an additional individual of the same species is 
created) with 50% probability provided the total energy of the 
system does not exceed Emax- When an empty energy level 
is picked, speciation occurs with 50% probability and a new 
individual of that species is created provided again the energy 
of the system does not exceed Emax- With 50% probability, 
no action is taken. This procedure is iterated until equilib- 
rium is reached. The effective temperature of the ecosystem 
is defined as the imbalance between Emax and the average 
energy of the system (Eq. Q). The figure shows a plot of 
the effective temperature of the ecosystem deduced from the 
simulations. The circles denote the data averaged over a run 
of 10^ time steps with the last 500 million used to compute 
the average temperature while the solid line is the theoretical 
prediction. 



We have carried out detailed computer simulations 
(Figs, n and [2Jl of systems with constraints on the to- 
tal energy and the total population and find excellent 
accord with theory. Fig. Q] shows the results of simula- 
tions corresponding to the case a ~ I with a constraint 
just on the total energy of the ecosystem. 

A strong hint that there is a link between the behav- 
iors of the ecosystem at Tmin and the physical system 
of Bosons at the BEC transition is obtained by noting 
that both these situations are characterized by a = 1 or 
equivalently a zero basic cost (/i = 0) for the introduction 
of an individual or a particle into the system. 

Let us set the value of N^ax equal to Nm , the average 
population in a system with Emax = Em and Nmax = oo, 
and consider the effect of varying the temperature. As in 



1.8 
1.6 
1.4 
1.2 
1 
0.8 
0.6 
0.4 
0.2 



1 1 1 r-1 1 1 1 — 






/ ^\. 












/ •*" 


.,.^9 








•"■ 


/ >"-"'"' 


.A. - 


_A 


/ * 






/ ^' 






f • 






/ • 






/ '' 






/ x 






/ A 






/ / 






• • 






/ ' 






/ -^ 












• / -^ 






/ ^ 






/ ^A 


















I 


. 


. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

T/T . 



FIG. 2: Phase transition in an ecosystem with Nmax ~ 
Nm = 65 and d = 1/2. The dashed and sohd curves are 
plots of theoretical predictions of jq-^ and ij, "'" respec- 
tively versus scaled temperature T /T^nin, where Tmm = 11-6 
and Em = 1000. The data points denote the results of simu- 



lations. d(E}/dT is a quantity analogous to the specific heat 
of a physical system and has the familiar A shape associated 
with the superfluid transition in liquid helium[y|. It was ob- 
tained in the simulations as the derivative of the interpolated 
values of (E) . The continuous phase transition is signaled by 
the peak in d{E)/dT (and the corresponding drop in (A'^)) 
on lowering the temperature and occurs in the vicinity of the 
temperature T,nin (the transition temperature moves closer 
to Train as the system size increases). 



Bose condensation, one can identify (for d > 0) a critical 
temperature, Tc, as the lowest temperature above which 



the first term on the rhs of Eq. H10() can be neglected[3- 
At Tc, {N) « Nmax = Nm and a « 1. Recall, however, 
that with only the energy ceiling, at T = Tmin the ecosys- 
tem was characterized by a = 1 and {N) — Nm- This 
confirms that Tc — Tmin- The macroscopic depletion of 
the community population, when T < Tmin is entirely 
akin to the macroscopic occupation of the ground state 
inBEC (Fig. El. 

When T is larger than Tmim the energy resources are 
sufficient to maintain the maximum allowed population 
and {N) = Nm- Physically, the state at Tmin corre- 
sponds to a maximally efficient use of the energy re- 
sources so that any decrease in T below Tmin inevitably 
leads to a decrease of the population and the biodiversity 
in the community. The existence of this novel transition 
is quite general and is independent of specific counting 
rules such as the ones used in the classic examples of 
distinguishable and indistinguishable particles Q. 

Our results suggest that ecosystems in equilibrium are 
organized in the vicinity of a novel critical point. Criti- 
cal points [13 are generally characterized by large fluctu- 
ations and a slowing down of dynamics in their vicinity, 
ft is an intriguing possibility that the ubiquity of power 
laws in ecology^, la, H Si ma-Y have some relationship to 
our findings here. 
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